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Abstract 

The leading diagrammatic two-loop corrections are incorporated into the prediction 
for the mass of the lightest Higgs boson, rrih, in the Minimal Supersymmetric Standard 
Model (MSSM). The results, containing the complete diagrammatic one- loop correc- 
tions, the new two-loop result and refinement terms incorporating leading electroweak 
two-loop and higher-order QCD contributions, are discussed and compared with results 
obtained by renormalization group calculations. Good agreement is found in the case 
of vanishing mixing in the scalar quark sector, while sizable deviations occur if squark 
mixing is taken into account. 



The search for the lightest Higgs boson provides a direct and very stringent test of 
Supersymmetry (SUSY), since the prediction of a relatively light Higgs boson is common 
to all Supersymmetric models whose couplings remain in the perturbative regime up to a 
very high energy scale IJ. A precise prediction for the mass of the lightest Higgs boson in 
terms of the relevant SUSY parameters is crucial in order to determine the discovery and 
exclusion potential of LEP2 and the upgraded Tevatron and also for physics at the LHC, 
where eventually a high-precision measurement of the mass of this particle might be possible. 

In the Minimal Supersymmetric Standard Model (MSSM) || the mass of the lightest 
Higgs boson, rrih, is restricted at the tree level to be smaller than the Z-boson mass. This 
bound, however, is strongly affected by the inclusion of radiative corrections. The dom- 
inant one-loop corrections arise form the top and scalar-top sector via terms of the form 
Gjrm! In (m^m^/mf) ||. They increase the predicted values of and yield an upper bound 
of about 150 GeV. These results have been improved by performing a complete one-loop 
calculation in the on-shell scheme, which takes into account the contributions of all sectors 
of the MSSM j|, ^, ||. Beyond one-loop order renormalization group (RG) methods have 
been applied in order to obtain leading logarithmic higher-order contributions []7|, |8], [], |TU| |, 



and a diagrammatic calculation of the dominant two-loop contributions in the limiting case 
of vanishing t-mixing and infinitely large Ma and tan/? has been carried out [IT]. 



The results of the latter calculations were found to yield considerably lower values for 
rrih than the one-loop on-shell calculation. The results obtained within the RG approach 
differ by up to 20 GeV from the one-loop on-shell result. The bulk of this difference can be 
attributed to the higher-order leading logarithmic contributions which are included in the 
RG results. The one-loop result on the other hand contains non-leading contributions which 
are not included in the RG results and whose size has not yet been precisely determined. 
Due to the large difference between the complete one-loop on-shell calculation and the RG 
results, and due to the difficulty in comparing the results of the two approaches, it is not easy 
to give an estimate for the accuracy of the current theoretical prediction for rrih in a similar 
way as done, for instance, for the electroweak precision observables within the Standard 
Model (SM) (see Ref. @). 

Recently a Feynman-diagrammatic calculation of the leading two-loop corrections of 
0(aa s ) to the masses of the neutral CP-even Higgs bosons has been performed |R|. Com- 



pared to the leading one-loop result the two-loop contribution was found to give rise to a 
considerable reduction of the rrih value. 

It is the purpose of this letter to combine the full diagrammatic one-loop on-shell result || 
with the leading diagrammatic two-loop result and to obtain in this way the currently 
most precise prediction for rrih within the Feynman-diagrammatic approach for arbitrary 
values of the parameters of the Higgs and scalar top sector of the MSSM. Further refinements 



concerning the leading two- loop Yukawa corrections of 0{G F rrq) || [14] and of leading 



QCD corrections beyond two-loop order are included. The resulting predictions for rrih are 
discussed in view of the discovery potential of LEP2 fl5| and are compared with the results 
obtained using the RG approach || ||, [T(J . 

Contrary to the SM, in the MSSM two Higgs doublets are needed. The Higgs potential 
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is given by |16| 

V = m 2 H l H l + m 2 2 H 2 H 2 -m 2 12 (e ab H 1 2 H 2 > + h.c.) 

+ {H X H X - H 2 H 2 f + ^ \H X H 2 \ 2 , (1) 

where m x ,m 2 ,m X2 are soft SUSY-breaking terms, g,g' are the SU(2) and £7(1) gauge cou- 
plings, and 612 = —1. 

The doublet fields H x and H 2 are decomposed in the following way: 

rr (Hi) _( V 1 + ( ( p° 1 +i X l )/V2\ 

H2 = (^) = («s + (^ + ix§)/^)- (2) 

The potential eq. ([I]) can be described with the help of two independent parameters (besides 
g, g')\ tan f3 = v 2 /v x and M\ = — m 2 2 (tan f3 + cot (3 ), where M^ is the mass of the CP-odd 
A boson. 

In order to obtain the CP-even neutral mass eigenstates, the rotation 



H° \ / cos a sin a \ ( (fi x 
h° J I — sin a cos a J I 4> 2 



(3) 



is performed, where the mixing angle a is given in terms of tan f3 and Ma by 

M\ + Ml 7T 

tan2a = tan2/3 M ^_ M f , -- < a < 0. (4) 

At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the (f) X — <p 2 
basis in terms of Mz and Ma by 



M 



2,tree = / 



Mi sin 2 /? + Mf cos 2 /3 — (M| + M|) sin/5 cos /? 
- {M\ + M| ) sin /3 cos (3 M\ cos 2 /3 + Mf sin 2 j3 



(5) 



which then has to be rotated with the angle a according to eq. (0), and one obtains the 
tree-level Higgs-boson masses 

Mffigss { m 2 jtre J- ^ 

In the Feynman-diagrammatic approach the one-loop corrected Higgs masses are derived 
by finding the poles of the h — //-propagator matrix whose inverse is given by 

/ A n-1 _ • /V -m 2 H:tvee + ± H (q 2 ) t hH (q 2 ) \ 

[ megs) { W? 2 ) q 2 -m 2 + ± h ( q 2 ) J' (7j 



where the £ denote the full one-loop contributions to the renormalized Higgs-boson self- 
energies. For these self-energies we take the result of the complete one-loop on-shell calcu- 
lation of Ref. IJ. The agreement with the result obtained in Ref. M is better than 1 GeV 
for almost the whole MSSM parameter space. 

As mentioned above the dominant contribution arises from the t — t-sector. The current 
eigenstates of the scalar quarks, and (Jr, mix to give the mass eigenstates q~i and q 2 . The 
non-diagonal entry in the scalar quark mass matrix is proportional to the mass of the quark 
and reads for the t-mass matrix 



m t M t LR = mt(A t - im cot P), 



where we have adopted the conventions used in Ref. [17]. Due to the large value of m t 
mixing effects have to be taken into account. Diagonalizing the t-mass matrix one obtains 
the eigenvalues mz and m~ t2 and the t mixing angle 9~ t . 

The leading two-loop corrections have been obtained in Ref. [13|| by calculating the 
O (ctctg) contribution of the t — t-sector to the renormalized Higgs-boson self-energies at 
zero external momentum from the Yukawa part of the theory. Since this calculation was 
performed in the <\>\ — 02-basis, we perform the rotation into the h — if -basis according to 
eq. ©: 



M2) 
p(2) 

-'h 

^(2) 
J hH 



2 v( 2 ) i • 2 v(2) . • ^(2) 

cos a + sin a L^ 2 + 2 sin a cos a 



• 2 v(2) 

sm a h\ 



2 v\(2) r, • v>( 2 ) 

cos a L^ 2 — 2 sm a cos a 2j^ 2 



sin a cos a 



^(2) 
J 4>2 



(cos 2 a 



■ 2 \v>(2) 

sm a)E^ 2 . 



(9) 



At the two-loop level the matrix (|7J) then consists of the renormalized Higgs-boson self- 
energies 

±.(f) = ±Wtf) + ±W{0), s = h,H,hH, (10) 

where the momentum dependence is neglected only in the two-loop contribution. The Higgs- 
boson masses at the two-loop level are obtained by determining the poles of the matrix AHi ggs 
in eq. flTJ). 

We have implemented two further steps of refinement into the prediction for mn, which are 
shown separately in the plots below. The leading two-loop Yukawa correction of 0(G 2 F mf) 
is taken over from the result obtained by renormalization group methods. It reads B O] 



Am 2 



with 
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The second step of refinement concerns leading QCD corrections beyond two-loop order, 
taken into account by using the MS top mass, m t = m t {m t ) f& 166.5 GeV, for the two-loop 
contributions instead of the pole mass, m t = 175 GeV. In the t mass matrix, however, we 
continue to use the pole mass as an input parameter. Only when performing the comparison 
with the RG results we use mt in the t mass matrix for the two- loop result, since in the RG 
results the running masses appear everywhere. This three-loop effect gives rise to a shift up 
to 1.5 GeV in the prediction for m^. 

For the numerical evaluation we have chosen two values for tan/3 which are favored 
by SUSY-GUT scenarios |TS[: tan/? = 1.6 for the SU(5) scenario and tan/3 = 40 for 



the 50(10) scenario. Other parameters are M z = 91.187 GeV, M w = 80.375 GeV, G F = 
1.16639 10 -5 GeV _2 ,a s (m t ) = 0.1095, and m t = 175 GeV. For the figures below we have 
furthermore chosen M = 400 GeV (M is the soft SUSY breaking parameter in the chargino 
and neutralino sector), Ma = 500 GeV, and m~ g = 500 GeV as typical values (if not indicated 
differently). The scalar top masses and the mixing angle are derived from the parameters 
M {l , M- tR and M^ R of the t mass matrix (our conventions are the same as in Ref. |I7| ). In 
the figures below we have chosen m g = M^ L = M% . 

The plot in Fig. [I] shows the result for m^ obtained from the diagrammatic calculation 
of the full one-loop and leading two-loop contributions. The two steps of refinement dis- 
cussed above are shown in separate curves. For comparison the pure one-loop result is also 
given. The results are plotted as a function of M^ R /m g , where m g is fixed to 500 GeV. 
The qualitative behavior is the same as in Ref. ||13|| , where the result containing only the 
leading one-loop contribution (and without further refinements) was shown. The two-loop 
contributions give rise to a large reduction of the one-loop result of 10-20 GeV. The two 
steps of refinement both increase mh by up to 2 GeV. A minimum occurs for M^ R = GeV 
which we refer to as 'no mixing'. A maximum in the two-loop result for is reached for 
about M^ R /m g « 2 in the tan/3 = 1.6 scenario as well as in the tan/3 = 40 scenario. This 
case we refer to as 'maximal mixing'. The maximum is shifted compared to its one-loop 
value of about M^ R jm g ~ 2.4. The two steps of refinement have only a negligible effect on 
the location of the maximum. 

In Fig. H] mh is shown in the two scenarios with tan/? = 1.6 and tan/3 = 40 as a 
function of m g for no mixing and maximal mixing. The tree-level, the one-loop and the 
two-loop results with the two steps of refinement are shown (the values of m g are such that 
the corresponding t-masses lie within the experimentally allowed region). In all scenarios of 
Fig. |2| the two-loop corrections give rise to a large reduction of the one-loop value of m^. The 
biggest effect occurs for the tan/? = 1.6 scenario with maximal mixing. The inclusion of 
the refinement terms leads to a slight shift in mh towards higher values, whose size is about 
20% of the two-loop correction. In the tan/3 = 1.6 scenario, mh reaches about 80 GeV for 
m g = 1 TeV in the no-mixing case, and about 100 GeV for m g = 1 TeV in the maximal- 
mixing case. For tan (3 = 40 the respective values of are nearly 115 GeV in the no- mixing 
case, and almost 130 GeV in the maximal-mixing case. 

Varying the gaugino parameter M, which enters via the non-leading one-loop contribu- 
tions, in the results of Fig. § changes the value of mh within 3 GeV. Different values of the 
gluino mass, m g , in the two-loop contribution affect the prediction for m h by up to 3 GeV. 
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Figure 1: One- and two-loop results for rrih as a function of M^ R /m,q for two values of tan /? 
The two steps of refinement discussed in the text are shown separately. 



Allowing for a splitting between the parameters M^ L , Mf in the t mass matrix yields max- 
imal values of rrth which are approximately the same as for the case rriq = M$ L = Mf (see 
also Fig. |5| below). Varying tan/? around the value tan/? = 1.6 leads to a relatively large 
effect in (higher values for rrih are obtained for larger tan /? ) , while the effect of varying 
tan p around tan /? = 40 is marginal. A more detailed analysis of the dependence of our 
results on the different SUSY parameters will be presented in a forthcoming publication. 
The discovery limit of LEP2 is expected to be slightly above 100 GeV fl5|]. Accordingly, 



our results confirm that for the scenario with tan/3 = 1.6 practically the whole parameter 
space of the MSSM can be covered at LEP2. For slightly larger tan/? and maximal mixing, 
however, some parameter space could remain in which the Higgs boson escapes detection at 
LEP2. For tan/? = 40, on the other hand, a full exploration of the MSSM parameter space 
will not be possible at LEP2. While the prediction for is at the edge of the LEP2 range 
in the no-mixing case, the case of large mixing will be a challenge for the upgraded Tevatron 
or will finally be probed at LHC. 

We now turn to the comparison of our diagrammatic results with the predictions obtained 
via renormalization group methods. We begin with the case of vanishing mixing in the 
t sector and large values of Ma, for which the RG approach is most easily applicable and is 
expected to work most accurately. In order to study different contributions separately, we 
have first compared the diagrammatic one-loop on-shell result || with the one-loop leading 
log result (without renormalization group improvement) given in Ref. [1C] and found very 
good agreement, typically within 1 GeV. We then performed a leading log expansion of our 
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Figure 2: The mass of the lightest Higgs boson for tan/3 = 1.6 and tan/3 = 40. The 
tree-level, the one-loop and the two-loop results for are shown as a function of rriq for 
the no-mixing and the maximal-mixing case. 



diagrammatic result (which corresponds to the two-loop contribution in the RG approach) 
and also found agreement with the full two-loop result within about 1 GeV. Finally, as 
shown in Fig. |3|, we have compared our diagrammatic result for the no-mixing case including 
the refinement terms with the RG results obtained in Ref. ||.[] As can be seen in Fig. [3], 
after the inclusion of the refinement terms the diagrammatic result for the no-mixing case 
agrees very well with the RG result. The deviation between the results exceeds 2 GeV only 
for tan/3 = 1.6 and trig < 150 GeV. For smaller values of Ma the comparison for the no- 
mixing case looks qualitatively the same. For tan/3 = 1.6 and values of Ma below 100 GeV 
slightly larger deviations are possible. Since the RG results do not contain the gluino mass 
as a parameter, varying rrig gives rise to an extra deviation, which in the no-mixing case 
does not exceed 1 GeV. Varying the other parameters /i and M in general does not lead to 
a sizable effect in the comparison with the corresponding RG results. 

We now consider the situation when mixing in the t sector is taken into account. We 
have again compared the full one-loop result with the one-loop leading log result used within 
the RG approach [|Kj and found good agreement. Only for values of Ma below 100 GeV and 
large mixing deviations of about 5 GeV occur. In Fig. |] our diagrammatic result including 
the refinement terms is compared with the RG results as a function of M^ R /mq for 
tan/? = 1.6 and tan/3 = 40. The point M^ R /rriq = corresponds to the plots shown in 
Fig. |3|, except that the parameter /i is now set to fi = —rriq. For larger t-mixing sizable 

1 The RG results of Ref. M and Ref. [Gil agree within about 2 GeV with each other. 
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Figure 3: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods |§. The mass of the lightest Higgs boson is 
shown for the two scenarios with tan/3 = 1.6 and tan/3 = 40 for the case of vanishing 
mixing in the t-sector. 



deviations between the diagrammatic and the RG results occur, which can exceed 5 GeV for 
moderate mixing and become very large for large values of M^ R /rrig. As already stressed 
above, the maximal value for rrih in the diagrammatic approach is reached for M^ R /rrig pa 2, 
whereas the RG results have a maximum at M^ R /rrig pa 2.4, i.e. at the one-loop value. 
Varying the value of rrig in our result leads to a larger effect than in the no-mixing case and 
shifts the diagrammatic result relative to the RG result within ±2 GeV. 

The results of our diagrammatic on-shell calculation and the RG methods have been 
compared above in terms of the parameters Mi , Mi and M^ R of the t mixing matrix, since 



the available numerical codes for the RG results || |10[ are given in terms of these parameters. 
However, since the two approaches rely on different renormalization schemes, the meaning of 
these (non-observable) parameters is not precisely the same in the two approaches starting 
from two-loop order. Indeed we have checked that assuming fixed values for the physical 
parameters , mi 2 , and &i and deriving the corresponding values of the parameters Mf , 
Mi R and M t LR in the on-shell scheme as well as in the MS scheme, sizable differences occur 
between the values of the mixing parameter M^ R in the two schemes, while the parameters 
Mi L , Mi are approximately equal in the two schemes. Thus, part of the different shape of 
the curves in Fig. |] may be attributed to a different meaning of the parameter M^ R in the 
on-shell scheme and in the RG calculation. 

In order to compare results obtained by different approaches making use of different 
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Figure 4: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods ||. The mass of the lightest Higgs boson is 
shown for the two scenarios with tan/3 = 1.6 and tan/3 = 40 for increasing mixing in the 
t-sector and rriq = Ma- 



renormalization schemes we find it preferable to compare predictions for physical observables 
in terms of other observables (instead of unphysical parameters). As a step into this direction 
we compare in Fig. [S] the diagrammatic results and the RG results as a function of the 
physical mass mf 2 and with the mass difference Am ( - = mf 2 — m>i x and the mixing angle 9~ t 
as parameters. In the context of the RG approach the running t masses, derived from the 
t mass matrix, are considered as an approximation for the physical masses. The range of 
the t masses appearing in Fig. |^ has been constrained by requiring that the contribution 
of the third generation of scalar quarks to the p-parameter |T7| does not exceed the value 
of 1.3 • 10~ 3 , which corresponds to the resolution of Ap (= ei) when it is determined from 
experimental data . As in the comparison performed above, in Fig. |5| very good agreement 



is found between the results of the two approaches in the case of vanishing t mixing. For 
the maximal mixing angle 6~ t = — 7r/4, however, the diagrammatic result yields values for 
which are higher by about 5 GeV. 

In summary, we have implemented the result of the Feynman-diagrammatic calculation 
of the leading 0(aa s ) corrections to the masses of the neutral CP-even Higgs bosons in the 
MSSM into the prediction based on the complete diagrammatic one-loop on-shell result. Two 
further refinements have been added in order to incorporate leading electroweak two-loop 
and higher-order QCD contributions. In this way we provide the at present most precise 
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Figure 5: Comparison between the Feynman-diagrammatic calculations and the results 
obtained by renormalization group methods ||. The mass of the lightest Higgs boson is 
shown for the two scenarios with tan/3 = 1.6 and tan/3 = 40 as a function of the heavier 
physical t mass m t - 2 . For the curves with 9f = a mass difference Amf = GeV is assumed 
whereas for 9$ = —tt/A we chose Am t - = 350 GeV, for which the maximal Higgs masses are 
achieved. 



prediction for rrih based on Feynman-diagrammatic results. The results are valid for arbitrary 
values of the relevant MSSM parameters. The two-loop corrections lead to a large reduction 
of the one-loop on-shell result. Concerning the discovery and exclusion potential of LEP2, 
our results confirm that the scenario with low tan/3 (tan/3 = 1.6) in the MSSM should be 
covered at LEP2, while the scenario with high tan (3 (tan (3 = 40) is only (partly) accessible 
for vanishing mixing in the scalar top sector. We have compared our results with the results 
obtained via RG methods, and have analyzed in this context the one-loop and two-loop 
contributions separately. We have found that the drastic deviations present between the 
one-loop on-shell result and the RG result are largely reduced. For the case of vanishing 
mixing in the scalar top sector the diagrammatic and the RG results agree very well, while 
sizable deviations exceeding 5 GeV occur when t mixing is taken into account. Since the 
gluino mass does not appear as a parameter in the RG results, its variation gives rise to an 
extra deviation which lies within ±2 GeV. We have discussed the issue of how the results 
obtained via different approaches using different renormalization schemes can be formulated 
such that they are readily comparable to each other also when corrections beyond one-loop 
order are incorporated. For this purpose it is very desirable to express the prediction for 
the Higgs-boson masses in terms of physical observables, i.e. the physical masses and mixing 
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angles of the model. 

We thank M. Carena, H. Haber and C. Wagner for fruitful discussions and communication 
about the numerical comparison of our results. We also thank A. Djouadi and M. Spira for 
helpful discussions. 
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